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The  existence  of  a  critical  roughness  height  for  outer  layer  similarity  between  smooth  and  rough 
wall  turbulent  boundary  layers  is  investigated.  Results  are  presented  for  boundary  layer 
measurements  on  flat  plates  covered  with  sandgrain  and  woven  mesh  with  the  ratio  of  the  boundary 
layer  thickness  to  roughness  height  ( S/k )  varying  from  16  to  110  at  Rew=7.3  X  103— 13  X  103.  In  all 
cases  tested,  the  layer  directly  modified  by  the  roughness  (the  roughness  sublayer)  is  confined  to  a 
region  <5 k  or  <3 ks  from  the  wall  (where  ks  is  the  equivalent  sandgrain  roughness  height).  In  the 
larger  roughness  cases,  this  region  of  turbulence  modification  extends  into  the  outer  flow.  However, 
beyond  5k  or  3 ks  from  the  wall,  similarity  in  the  turbulence  quantities  is  observed  between  the 
smooth  and  rough  wall  boundary  layers.  These  results  indicate  that  a  critical  roughness  height, 
where  the  roughness  begins  to  affect  most  or  all  of  the  boundary  layer,  does  not  exist.  Instead,  the 
outer  flow  is  only  gradually  modified  with  increasing  roughness  height  as  the  roughness  sublayer 
begins  to  occupy  an  ever  increasing  fraction  of  the  outer  layer.  [DOI:  10.1063/1.2757708] 


I.  INTRODUCTION 

Smooth  wall  studies  have  provided  the  basis  for  under¬ 
standing  rough  wall  flows,  with  engineering  models  gener¬ 
ally  treating  surface  roughness  as  a  small  perturbation  to  the 
smooth  wall  boundary  layer.  However,  questions  remain  as 
to  the  extent  of  the  roughness-induced  perturbation  and 
whether  or  not  smooth  wall  boundary  layer  theory  is  valid 
for  large  roughness.  A  fundamental  question  with  implica¬ 
tions  for  the  use  of  wall  models  and  large  eddy  simulations 
for  the  computation  of  rough  wall  flows  is  if  the  buffer  layer 
can  be  perturbed  or  destroyed  without  transmitting  anything 
to  the  outer  flow  except  a  change  in  skin  friction. 1  If  so,  what 
happens  for  larger  roughness  elements? 

The  major  effect  of  roughness  on  the  mean  velocity  pro¬ 
file  plotted  in  inner  variables  is  a  downward  shift  in  the  log- 
law  region  (A  U+),  termed  the  roughness  function,  while  the 
shape  of  the  profile  in  the  outer  layer  is  unchanged.  In  outer 
scaling,  this  results  in  the  collapse  of  smooth  and  rough  pro¬ 
files  in  velocity-defect  form,  which  is  supported  by  numer- 
ous  studies.  This  mean  flow  similarity  indicates  that  the 
use  of  wall  functions  for  rough  wall  boundary  layer  models 
may  be  appropriate.  However,  results  from  Krogstad  et  al.9 
and  Keirsbulck  et  al. 10  indicate  that  the  wake  strength  of  the 
mean  velocity  profile  is  increased  on  rough  walls,  implying 
the  effect  of  surface  roughness  may  propagate  well  into  the 
outer  region.  If  this  is  the  case,  application  of  rough  wall 
models  that  rely  on  outer  layer  similarity  will  yield  erroneous 
results. 

Townsend’s11  Reynolds  number  similarity  hypothesis 
and  subsequent  extensions  by  Perry  and  Chong  and 
Raupach  et  al.4  state  that  the  turbulent  stresses  are  indepen¬ 
dent  of  surface  condition  outside  the  roughness  sublayer  (or 
viscous  sublayer  in  the  smooth  wall  case)  at  sufficiently  high 


Reynolds  number  when  normalized  by  the  wall  shear  stress. 
That  is,  the  turbulence  in  the  outer  layer  is  unaffected  by 
surface  roughness  except  in  the  role  it  plays  in  modifying  the 
outer  velocity  and  length  scales.  The  underlying  assumption 
of  Townsend’s  hypothesis  is  that  the  boundary  layer  thick¬ 
ness  (5)  is  large  compared  to  the  roughness  height  (k).  The 
roughness  sublayer  is  typically  thought  to  extend  ~5  rough¬ 
ness  heights  from  the  wall.4  A  number  of  studies  have  shown 
support  for  Townsend’s  hypothesis  for  three-dimensional 
roughness  including  the  works  of  Ligrani  and  Moffat  for 
packed  spheres,  Perry  and  Li14  for  expanded  mesh,  Kunkel 
and  Marusic15  for  salt  flats  ground  cover,  Shockling  et  al . 16 
for  a  honed  pipe,  and  the  present  authors  for  sanded 
surfaces,6  packed  spheres  with  and  without  grit,7  and  sand- 
paper  and  woven  mesh.  Recently,  the  present  authors  also 
demonstrated  support  for  the  hypothesis  for  the  limiting  case 
of  fully  rough  flow  at  high  Reynolds  number  for  a  roughness 
that  is  a  very  small  fraction  of  the  boundary  layer  thickness. 

In  contrast,  changes  to  the  Reynolds  stresses  over  a  large 
portion  of  the  outer  layer  have  been  observed  for  flows  over 
woven  mesh  and  transverse  bar  roughness  in  the  experimen¬ 
tal  work  of  Krogstad  and  Antonia,18  Tachie  et  al.,19  and 
Keirsbulck  et  al. 10  Recent  numerical  simulations  of  turbulent 
channel  flow"  ’  have  also  shown  roughness  effects  extend¬ 
ing  into  the  outer  layer.  A  possible  explanation  for  the  dis¬ 
parate  findings  may  be  due  to  the  “strong”  roughness  used  in 
the  investigations,  where  effects  were  observed  in  the  outer 
layer.  A  strong  roughness  can  be  thought  of  as  a  surface 
whose  roughness  height  (or  alternatively  equivalent 
sandgrain  roughness  height")  is  a  significant  portion  of  the 
inner  layer  thickness,  but  just  how  significant  a  portion  is 
unclear.  Therefore,  the  maximum  roughness  height  for  which 
outer  layer  similarity  may  be  assumed  to  be  valid  needs  to  be 
defined.  For  cases  with  larger  roughness  than  this,  the  clas- 
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sical  concept  of  the  roughness  being  a  small  perturbation  to 
the  smooth  wall  boundary  layer,  with  its  effect  characterized 
solely  by  AU+,  must  be  abandoned.  Based  upon  his  review, 
Jimenez1  suggested  that  similarity  in  the  outer  flow  can  be 
expected  for  rough  wall  flows  with  S/k  as  a  minimum 
greater  than  40,  but  likely  the  limiting  value  is  closer  to  80. 
Recently,  Schultz  and  Flack  and  Flack  et  al.  have  observed 
similarity  for  roughness  with  S/k  of  29  and  45,  respectively. 
The  proper  roughness  length  scale  will  likely  be  more  than 
simply  a  measure  of  the  roughness  height  to  boundary  layer 
thickness,  but  moreover  will  include  the  effects  of  the  rough¬ 
ness  type  and  density  on  the  mean  and  turbulent  flow.  Using 
the  equivalent  sandgrain  roughness  as  a  roughness  scale  does 
provide  a  common  currency  for  a  variety  of  roughness 
types,  since  it  incorporates  the  effect  of  the  surface  rough¬ 
ness  on  the  mean  flow.  However,  it  may  not  be  the  best 
indicator  of  effects  on  the  Reynolds  stresses  and  higher  order 
moments,  and  has  the  disadvantage  of  not  being  known  for  a 
generic  roughness  a  priori. 

The  goal  of  the  present  research  is  to  document  the  ef¬ 
fect  of  increasing  roughness  height  on  the  outer  layer  turbu¬ 
lence  statistics  in  fully  developed  turbulent  boundary  layers 
on  two  types  of  three-dimensional  roughness.  The  range  of 
S/k  was  16=5  S/k^  110  (6^  S/ ks^9l)  in  the  study.  This 
range  spans  from  where  outer  layer  similarity  is  expected  to 
be  valid  to  significantly  rougher  surfaces,  where  similarity  is 
expected  to  fail  based  on  the  criteria  of  Jimenez. 1  From  these 
results,  the  similarity  of  both  the  mean  velocity  and  turbu¬ 
lence  structure  on  rough  and  smooth  walls  is  critically  evalu¬ 
ated.  Additionally,  the  limiting  roughness  height  where  outer 
layer  similarity  holds  is  explored. 


II.  EXPERIMENTAL  FACILITIES  AND  METHODS 


The  experiments  were  conducted  in  the  U.S.  Naval 
Academy’s  large  recirculating  water  tunnel.  The  test  section 
is  40  X  40  cm2  in  cross  section  and  is  1.8  m  in  length,  with  a 
tunnel  velocity  range  of  0-8.0  m/s  and  a  freestream  turbu¬ 
lence  intensity  of  —0.5%.  The  present  tests  were  run  at 
a  tunnel  speed  of  —2.0  m/s  (Rev=2.7  X 106,  Refl=6.1 
X  103-13X  103).  The  test  plates  were  flush  mounted  into  a 
splitter-plate  test  fixture.  The  test  fixture  was  mounted  hori¬ 
zontally,  at  mid-depth  in  the  tunnel.  The  first  200  mm  of  the 
test  fixture  is  covered  with  #36-grit  sandpaper  to  ensure  ad¬ 
equate  turbulent  boundary  layer  tripping.  Measurements 
were  obtained  1.35  m  downstream  of  the  trip,  allowing  for  a 
sufficient  boundary  layer  growth.  The  upper,  removable  wall 
of  the  tunnel  is  adjustable  to  account  for  boundary  layer 
growth,  producing  a  nearly  zero  pressure  gradient  boundary 
layer,  with  an  acceleration  parameter  of  1  X  10"8,  where 
K  is  defined  as 


v  dUe 
U2e  dx 


(1) 


A  plan  view  schematic  of  the  facility  is  shown  in  Fig.  1. 
Additional  details  of  the  experimental  facility  can  be  found 
in  Schultz  and  Flack6  and  Flack  et  al* 

Seven  surfaces  were  tested  in  this  study,  as  listed  in 
Table  I:  one  was  a  smooth  cast  acrylic  surface;  the  other  six 


sandpaper  trip  test  surface  measurement  location 


FIG.  1.  Schematic  of  flat  plate  test  fixture. 


were  rough  surfaces.  Three  surfaces  were  covered  with  wet/ 
dry  sandpaper.  These  included  80-,  24-,  and  12-grit  sandpa¬ 
per.  The  remaining  three  were  covered  with  woven  wire 
mesh  with  pitch  to  diameter  ratios  of  6.25  for  Ml,  4.58  for 
M2,  and  8.45  for  M3.  The  roughness  height,  k,  given  in 
Table  1  is  the  maximum  peak  to  trough  height  measured  over 
a  50  X  50  mm2  sampling  area.  For  the  sandpaper  surfaces, 
this  was  measured  using  a  dial  indicator.  For  the  mesh  sur¬ 
faces,  it  was  taken  to  be  twice  the  wire  diameter.  The  correct 
scale  to  use  to  accurately  characterize  the  roughness  height  is 
still  an  active  area  of  research  and  may  depend  on  additional 
parameters  to  account  for  differences  in  texture. 

Velocity  measurements  were  made  using  a  TSI  FSA3500 
two-component,  fiber-optic  laser  Doppler  velocimeter 
(LDV).  The  LDV  consists  of  a  four-beam  fiber-optic  probe 
that  collects  data  in  backscatter  mode.  A  custom  designed 
beam  displacer  was  added  to  the  probe  to  shift  one  of  the 
four  beams,  resulting  in  three  coplanar  beams  that  can  be 
aligned  parallel  to  the  wall.  This  allowed  for  near-wall  mea¬ 
surements  without  having  to  tilt  the  probe  at  a  small  angle  or 
rotate  the  probe  to  resolve  velocity  components.  The  mea¬ 
surement  volume  diameter  was  90  pm  and  its  length  was 
1.3  mm.  The  viscous  length  scale  (u/«T)  varied  from  7.5  pm 
for  the  roughest  surface  to  13  pm  for  the  smooth  wall.  The 
diameter  of  the  probe  volume,  therefore,  ranged  from  7  to  1 3 
viscous  lengths  in  the  present  study.  The  flow  was  seeded 
with  2  pm  alumina  particles.  The  seed  volume  was  con¬ 
trolled  to  achieve  acceptable  data  rates  while  maintaining  a 
low  burst  density  signal."  A  total  of  40  000  random  velocity 
samples  were  obtained  in  coincidence  mode  at  each  location 
in  the  boundary  layer.  Doppler  bursts  for  the  two  channels 
were  required  to  fall  within  a  50- ps  coincidence  window  or 
the  sample  was  rejected.  The  probe  was  traversed  to  approxi¬ 
mately  40  locations  within  the  boundary  layer  with  a  Velmex 
three-axis  traverse  unit.  The  traverse  allowed  the  position  of 
the  probe  to  be  maintained  to  +5  pm  in  all  directions. 

Precision  uncertainty  estimates  for  the  velocity  measure¬ 
ments  were  made  through  repeatability  tests  using  the  proce- 
dure  given  by  Moffat.  Ten  replicate  velocity  profiles  were 
taken  on  both  a  smooth  and  a  rough  plate.  LDV  measure¬ 
ments  are  susceptible  to  a  variety  of  bias  errors  including 
angle,  velocity,  and  velocity  gradient  biases,  as  detailed  by 
Edwards."6  Fringe  bias  results  from  the  inability  to  sample 
scattering  particles  passing  through  the  measurement  volume 
at  large  angles  since  several  fringe  crossings  are  needed  to 
validate  a  measurement.  In  this  experiment,  the  fringe  bias 
was  considered  insignificant,  as  the  beams  were  shifted  well 
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TABLE  I.  Experimental  test  conditions. 


Surface 

Symbol 

u„ 

(m/s) 

Re„ 

Wall  offset 
(mm) 

uT  Clauser 
(m/s) 

utTS 

(m/s) 

A  U* 

8 

(mm) 

<5* 

(mm) 

e 

(mm) 

k 

(mm) 

k+ 

K 

S/k 

Slks 

Smooth 

2.00 

6140 

0.0763 

0.0746 

32.5 

2480 

4.05 

3.08 

80-grit 

sandpaper 

SI 

1.98 

7970 

0.175 

0.0917 

0.0901 

5.24 

37.2 

3250 

5.93 

4.22 

0.690 

60 

36 

54 

91 

Fine 

mesh 

Ml 

1.99 

7290 

0.130 

0.0981 

0.0951 

6.31 

34.8 

3080 

5.83 

4.06 

0.320 

28 

56 

110 

55 

24-grit 

sandpaper 

S2 

1.99 

8970 

0.580 

0.120 

0.119 

10.8 

37.9 

4140 

7.91 

4.96 

1.80 

200 

360 

21 

11 

Medium 

mesh 

M2 

2.00 

9110 

0.400 

0.119 

0.113 

10.9 

38.3 

4180 

8.00 

5.00 

1.40 

150 

370 

27 

11 

12-grit 

sandpaper 

S3 

2.00 

12  430 

1.130 

0.132 

0.131 

13.0 

44.9 

6060 

9.99 

6.08 

2.85 

380 

860 

16 

7.0 

Coarse 

mesh 

M3 

2.00 

13  050 

0.810 

0.133 

0.126 

13.7 

46.8 

6350 

10.91 

6.40 

2.45 

330 

1150 

19 

5.5 

above  a  burst  frequency  representative  of  twice  the 
freestream  velocity.  Validation  bias  results  from  filtering  too 
near  the  signal  frequency  and  any  processor  biases.  In  gen¬ 
eral  these  errors  are  difficult  to  estimate  and  vary  from  sys¬ 
tem  to  system.  No  corrections  were  made  to  account  for 
validation  bias.  Velocity  bias  results  from  the  greater  likeli¬ 
hood  of  high  velocity  particles  moving  through  the  measure¬ 
ment  volume  during  a  given  sampling  period.  The  present 
measurements  were  burst  transit  time  weighted  to  correct  for 
velocity  bias,  as  presented  in  Buchhave  et  al.  Velocity  gra¬ 
dient  bias  is  due  to  velocity  variation  across  the  measurement 
volume.  These  bias  errors  only  become  significant  very  near 
the  wall,  and  in  the  present  study  these  were  quite  small,  and 
therefore,  no  velocity  gradient  bias  corrections  were  made. 
Bias  estimates  were  combined  with  the  precision  uncertain¬ 
ties  to  calculate  the  overall  uncertainties  for  the  measured 
quantities,  as  listed  with  the  figures  of  the  reduced  data. 

Two  methods  were  used  to  determine  the  friction  veloc¬ 
ity,  for  both  the  smooth  and  rough  surfaces,  as  listed  in 
Table  I.  The  first  was  the  total  stress  method,  which  assumes 
that  a  nominally  constant  shear  stress  region  exists  in  the 
inner  part  of  the  boundary  layer  which  is  equal  to  the  wall 
shear  stress.  The  total  stress  was  calculated  at  the  plateau  of 
the  Reynolds  shear  stress  profile  in  the  overlap  region  of  the 
boundary  layer  by  summing  the  contributions  of  the  viscous 
and  turbulent  stresses.  For  the  smooth  wall,  the  friction  ve- 
locity  was  also  determined  using  the  Clauser”  chart  method, 
with  log-law  constants  k=0.41  and  B  =  5.0.  A  modified 
Clauser  chart  method  was  also  employed  to  determine  the 
friction  velocity  on  the  rough  wall.14'17  This  method  also 
solves  for  the  error  in  origin,  or  wall  offset  due  to  roughness, 
as  listed  in  Table  I.  The  friction  velocity  values  obtained 
using  the  Clauser  (smooth)  and  modified  Clauser  (rough) 
methods  are  used  in  the  subsequent  data  reduction;  however, 
using  uT  determined  from  the  total  stress  method  would  not 
change  the  trends  in  the  presented  results  since  the  values 
obtained  from  both  methods  show  good  agreement.  The  un¬ 
certainty  in  uT  for  the  smooth  walls  using  the  Clauser  chart 
method  is  ±3%,  and  the  uncertainty  in  uT  for  the  rough  walls 
using  the  modified  Clauser  chart  method  was  +5%.  The  un¬ 


certainty  in  uT  using  the  total  stress  method  is  +7%  for  both 
the  smooth  and  rough  walls. 

The  boundary  layer  parameters  S,  cT,  S\  and  0  for  the 
smooth  and  rough  walls  are  listed  in  Table  I.  The  physical 
growth  of  the  boundary  layer  was  significantly  higher  on  the 
rough  walls,  with  the  displacement  and  momentum  thick¬ 
nesses  more  than  doubling  for  the  roughest  surface  as  com¬ 
pared  to  the  smooth  wall  for  the  same  freestream  velocity. 
For  this  study,  S  is  taken  as  the  wall  normal  location  where 
the  mean  velocity  is  99.5%  of  the  freestream  velocity.  The 
values  <5*  and  6  were  calculated  assuming  a  linear  velocity 
profile  from  the  data  point  closest  to  the  wall  and  v  =  0  for  the 
smooth  surface.  For  the  rough  surfaces,  a  logarithmic  veloc¬ 
ity  profile  from  the  data  point  closest  to  the  wall  and  y  =  0 
was  assumed. 

III.  RESULTS  AND  DISCUSSION 
A.  Mean  flow 

Mean  velocity  profiles  for  the  seven  surfaces  plotted  in 
inner  variables  are  shown  in  Fig.  2.  The  smooth  wall  results 
of  DeGraaff  and  Eaton”  taken  at  Re w= 5200  are  shown  for 
comparison.  Both  types  of  rough  surfaces  display  a  linear- 


FIG.  2.  Mean  velocity  profiles  for  all  surfaces  in  inner  variables.  U*± 3% 
for  smooth  wall  and  +5%  for  rough  walls. 
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FIG.  3.  Mean  velocity  profiles  in  velocity  defect  form  for  all  surfaces:  (a) 
classic  scaling;  (b)  Rotta-Clauser  scaling.  (Ue—U)IUT±3%  for  smooth  wall 
and  +5%  for  rough  walls. 

log  region  that  is  shifted  by  the  roughness  function,  Af/+, 
below  the  smooth  profile,  indicating  an  increased  momentum 
deficit  on  these  surfaces.  The  sandgrain  and  mesh  surface 
produced  roughness  functions  range  from  Af/+=5.2  to  13.7, 
as  listed  in  Table  I.  It  is  of  note  that  similar  roughness  func¬ 
tions  are  obtained  for  the  mesh  and  sandgrain  surfaces  at 
nominally  the  same  unit  Reynolds  number  even  though  the 
roughness  heights,  k,  differ  substantially.  This  is  most  evi¬ 
dent  for  the  80-grit  sandpaper  and  fine  mesh.  The  mesh  pro¬ 
duces  a  roughness  function  20%  larger  than  the  sandpaper 
although  the  roughness  height  of  the  sandpaper  is  more  than 
twice  that  of  the  mesh.  The  equivalent  sand  roughness 
height,  ks ,  defined  as  the  uniform  sand  roughness  height  that 
gives  that  same  roughness  function  in  the  fully  rough 
regime,  is  also  listed  in  Table  I.  ks  provides  a  common 
measure  of  the  influence  of  the  roughness  on  the  mean  flow 
and  is  related  to  Af/+  as  shown  in  the  following: 

AU+=-\n(k*)  +  B-C,  (2) 

K 

where  C  is  the  roughness  function  intercept  for  uniform  sand 
roughness  equal  to  8.5. 

The  mean  velocity  profiles  for  all  test  cases  are  pre¬ 
sented  in  velocity-defect  form  in  Fig.  3.  Again,  shown  for 
comparison  are  the  smooth  wall  results  of  DeGraaff  and 
Eaton.  '0  Figure  3(a)  shows  the  results  using  classic  outer 
scaling  with  normalization  by  uT  and  S.  Excellent  collapse  of 
the  present  smooth  and  rough  wall  results  is  seen  in  the  outer 


part  of  the  boundary  layer.  The  collapse  is  slightly  improved 
using  the  Rotta-Clauser  length  scale  ( <f  U+),  as  shown  in  Fig. 
3(b).  These  results  indicate  that  the  outer  layer  mean  flow  is 
largely  insensitive  to  surface  condition.  Further  discussion  of 
the  mean  flow  profiles  and  velocity-defect  scaling  for  the 

3 1 

present  surfaces  is  found  in  Connelly  et  al. 

B.  Reynolds  stresses  and  quadrant  analysis 

The  Reynolds  stress  profiles  for  the  rough  and  smooth 
surfaces  will  now  be  presented.  Ideally,  the  Karman  number 
(cT)  should  be  matched  to  investigate  outer  layer  similarity 
between  smooth  and  rough  surfaces.  If  the  disparity  between 
the  smooth  and  rough  wall  Karman  number  is  too  large  it 
would  be  difficult  to  discern  if  changes  in  the  profiles  in  the 
outer  layer  are  due  to  roughness  or  Reynolds  number  effects. 
It  has  been  shown  that  the  profiles  of  streamwise  Reynolds 
normal  stress,  u'2  ,  exhibit  an  increasing  peak  in  the  buffer 
layer  (at  y+~  15)  with  increasing  Karman  number  for  smooth 
wall  flows.  A  more  slowly  rising  log-layer  plateau  in  u'2 

29 

has  also  been  noted  in  the  studies  of  DeGraaff  and  Eaton' 

32 

and  Fernholz  and  Finley  for  boundary  layers,  and  Morrison 
et  al:  for  fully  developed  pipe  flow.  Figure  4  shows  smooth 
wall  Reynolds  stress  profiles  for  a  range  of  ^  similar  to  that 
covered  in  this  study.  Included  in  this  graph  are  the  smooth 
wall  results  from  the  current  study  (/T=2480),  a  higher 
Karman  number  smooth  wall  case  obtained  in  the  same 
facility17  (<5+  =  5170),  and  the  smooth  wall  results  of  DeGraff 
and  Eaton29  ( rT  =  1780  and  4570).  It  should  be  noted  that  the 
apparent  differences  observed  in  all  the  Reynolds  stresses  in 
the  near-wall  region  ( y  /  <5<0.03)  are  primarily  due  to  the  use 
of  outer  scaling.  There  is  a  slight  increase  in  the  log-layer 
plateau  for  the  streamwise  normal  Reynolds  stress  with  in¬ 
creasing  Karman  number  [Fig.  4(a)]  at  0.03  <y/ <-? < 0 . 1 .  Sig¬ 
nificant  differences  are  not  observed  in  u'2  for  \7 <5>0. 1 . 
Reasonable  agreement  is  observed  for  v' 2  and  -u'u'  in  the 
overlap  and  outer  layer.  There  do  not  appear  to  be  any  de¬ 
finitive  trends  in  these  Reynolds  stresses  with  increasing 
Karman  number  in  agreement  with  the  conclusions  of  De- 
Graaff  and  Eaton."  The  Karman-number-dependent  differ¬ 
ences  in  the  Reynolds  stresses  in  the  overlap  and  outer  flow 
are  only  observed  in  u'2  and  are  small  (well  within  the 
experimental  uncertainty).  Therefore,  subsequent  compari¬ 
sons  of  the  rough  wall  profiles  will  be  made  to  the  present 
smooth  wall  profiles  at  (5+=2480. 

Figures  5-7  show  the  Reynolds+normalized  Reynolds 
stress  profiles  ( u '2  ,  v'2  ,  and  -u'v'  )  for  the  smooth  and 
rough  walls.  In  the  first  graph  (a),  all  the  surfaces  are  pre¬ 
sented  together,  along  with  the  smooth  wall  results  of  De¬ 
Graaff  and  Eaton,29  at  Rctt=5200  ( <f  =  1 780),  for  compari¬ 
son.  Additionally,  smooth  and  rough  wall  results  are  shown 
in  log  scaling  to  highlight  the  approximate  location  where 
the  profiles  diverge.  The  lines  on  Figs.  5-7  represent  5k/  S 
(solid  line)  and  3 kj  8  (dashed  line)  for  each  of  the  rough 
surfaces,  with  endpoint  symbols  identifying  the  type  of  sur¬ 
face.  Five  roughness  heights4  and/or  three  equivalent 
sandgrain  roughness  heights8  have  previously  been  identified 
as  potential  wall-normal  limits  of  the  roughness  sublayer. 
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FIG.  4.  Smooth  wall  Reynolds  stress  profiles  in  outer  variables  for  a  range 
of  (a)  u’2  ;  (b)  v'2  ;  (c)  —u'u'  .  Symbols  listed  in  (a). 


The  streamwise  Reynolds  normal  stresses,  u'2  (Fig.  5) 
for  the  smooth  and  rough  walls  show  excellent  agreement 
throughout  the  outer +region  of  the  boundary  layer.  Outer 
layer  similarity  of  u'2  has  been  observed  in  numerous  other 
studies  of  rough  wall  boundary  layers,  as  discussed  in  the 
review  articles  of  Raupach  et  al ,4  and  Jimenez.1  Near  the 
wall  {yl  8^  0.04),  u ' 2  is  significantly  lower  on  the  rough 
walls.  For  the  mildest  roughness  (Ml,  SI),  this  corresponds 
to  y+;=  120-13,  near  the  location  of  the  outer  edge  of  the 
buffer  layer  in  a  smooth  wall  flow.  The  suppression  of  u'2 
increases  with  increasing  roughness  Reynolds  number  for 
flows  in  the  transitionally  rough  regime.  For  flows  in  the 
fully  rough  regime,  the  near-wall  peak  in  u'2  is  absent  en- 
tirely.  Ligrani  and  Moffat  ~  identified  the  absence  of  a  near¬ 
wall  peak  in  u'2  as  an  indicator  of  fully  rough  flow.  Flores 
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FIG.  5.  Streamwise  Reynolds  normal  stress  profiles  in  outer  variables:  (a) 
all  profiles;  (b)  smooth,  Ml  and  SI:  (c)  smooth.  M2,  and  S2;  (d)  smooth, 

M3,  and  S3.  ( — )  5k / S;  ( — )  3 ks/S;  uncertainty  in  u'2  +8%.  Symbols  listed 
in  (a). 


and  Jimenez34  elucidated  the  physical  reasons  for  this  obser¬ 
vation  identifying  that  the  low-speed  streaks  and  quasi- 
streamwise  vortices  of  the  buffer  layer  are  shortened  in 
rough  wall  flows.  This,  in  turn,  results  in  decreased  stream- 
wise  velocity  fluctuations.  For  fully  rough  flows,  surface 
roughness  eliminates  near-wall  viscous  effects,  necessarily 
altering  the  entire  viscous  sublayer.  Therefore,  the  likely 
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FIG.  6.  Wall-normal  Reynolds  normal  stress  profiles  in  outer  variables:  (a) 
all  profiles;  (b)  smooth,  Ml,  and  SI;  (c)  smooth,  M2,  and  S2;  (d)  smooth, 

M3,  and  S3.  ( — )  5k/ <5;  ( — )  3kJ r);  uncertainty  in  v'2  +9%.  Symbols  listed 
in  (a). 

minimum  wall-normal  extent  of  roughness  influence  for  fully 
rough  flows  is  the  location  of  the  inner  limit  of  the  log-law 
region,  even  for  small  roughness  heights. 

_ + 

The  wall-normal  Reynolds  normal  stresses,  v'2  (Fig.  6) 
also  indicate  good  collapse  in  the  outer  region  of  the  bound¬ 
ary  layer.  This  result  is  in  agreement  with  a  number  of 
studies.6-8'13,14  In  contrast,  the  studies  of  Krogstad  et  al.9  and 


y/8 


FIG.  7.  Reynolds  shear  stress  profiles  in  outer  variables:  (a)  all  profiles;  (b) 
smooth,  Ml,  and  SI;  (c)  smooth,  M2,  and  S2;  (d)  smooth,  M3,  and  S3.  ( — ) 

5 k/ S',  ( — )  3 ks! S\  uncertainty  in  -u'v'  ±9%.  Symbols  listed  in  (a). 


Keirsbulck  et  al.10  showed  that  significant  increases  in  the 
wall-normal  Reynolds  normal  stress  penetrate  well  into  the 
outer  layer  over  rough  walls.  These  findings  indicated  that 
there  may  be  a  critical  roughness  height  where  wall  similar¬ 
ity  fails,  prompting  the  systematic  increase  of  roughness 
height  in  the  present  study.  The  roughness  in  the  Krogstad 
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et  al.9  and  Keirsbulck  et  al.10  studies  had  S/ks=  15  and  7, 
respectively;  values  similar  to  the  larger  mesh  and  sandpaper 
roughness  of  the  present  study. 

_ + 

The  normalized  Reynolds  shear  stress  (- u'v '  )  profiles 
are  presented  in  Fig.  7.  Good  collapse  of  the  rough  and 
smooth  profiles  is  again  observed  in  the  outer  region  of  the 
boundary  layer  outside  of  a  roughness  sublayer.  While  near¬ 
wall  roughness  effects  are  present  for  all  the  Reynolds 
stresses,  similarity  exists  between  rough  and  smooth  walls  in 
the  outer  layer,  outside  of  the+roughness  sublayer.  The  agree¬ 
ment  of  the  v'2  and  —u'v'  profiles  supports  the  concept 
of  universal  active  motions  in  the  outer  layer,  where  the  ac¬ 
tive  motions  are  the  turbulent  motions  that  contribute  to  the 
Reynolds  shear  stress. 

The  present  results  are  consistent  with  previous  research 
that  has  shown  similarity  outside  a  roughness  sublayer  ex¬ 
tending  5k  or  3 ks  from  the  wall.  As  observed  in  Figs.  5-7, 
the  Reynolds  stresses  show  good  collapse  with  smooth  wall 
results  for  y^5k  or  y^3ks.  For  the  roughest  surfaces  (M3, 
S3),  5k  and  3 ks  seem  to  overpredict  the  region  of  roughness 
influence,  with  a  better  estimate  being  2k  or  1  ks.  For  the 
moderate  (M2,  S2)  and  mild  roughness  (Ml,  SI),  5k  and  3 ks 
appear  to  be  appropriate. 

Changes  in  turbulence  structure  between  rough  and 
smooth  walls  are  investigated  using  standard  quadrant  de- 
composition  techniques,  as  described  in  Wallace  et  air  us¬ 
ing  hyperbolic  hole  size,36  H,  in  the  trigger  function.  The 
contribution  to  u'v'  from  a  given  quadrant,  Q,  can  be  ex¬ 
pressed  as 


—  i  r 

( u'v  r  )Q=  lim  -  u'v' (t)IQ{t)dt,  (3) 

1  J  o 

where  Ig(t)  is  a  trigger  function  defined  as 

[l  when  \u'v'\o  2=  H(u'2)l,2{v'2)112 
Iq=\  Q  .  (4) 

[0  otherwise. 

Figure  8  shows  the  contributions  to  the  conditionally  aver¬ 
aged  Reynolds  shear  stress  for  all  surfaces  for  the  case  with 
a  hyperbolic  hole  size  H= 0  and  identifying  stronger  events 
with  H=  2.  For  H=  0  there  is  excellent  collapse  in  the  outer 
layer  between  the  rough  and  smooth  surfaces.  When  stronger 
events  are  identified  using  H= 2,  again  the  similarity  between 
rough  and  smooth  walls  is  within  the  uncertainty  of  the  mea¬ 
surement,  with  more  near-wall  scatter  for  the  ejection  events 
(Q2).  Krogstad  et  al.9  observed  greater  contributions  to  the 
Reynolds  shear  stress  from  the  Q2  and  Q4  quadrants  for  a 
rough  wall  covered  with  a  similar  mesh  surface  ( S/ks - 15)  as 
compared  to  a  smooth  wall.  Outer  layer  collapse  was  ob¬ 
served  for  the  Q4  events  when  the  strongest  sweeps  were 
identified  with  a  hyperbolic  hole  size  of  H= 4.  Schultz  and 
Flack7  observed  stronger  Q4  and  Q2  events  near  the  wall  for 
flow  over  a  packed  sphere  bed,  due  to  the  less  strict  bound¬ 
ary  condition  for  a  rough  wall.  The  contributions  from  the 
ejection  and  sweep  events  collapsed  with  the  smooth  wall 
profiles  in  the  outer  layer.  Similar  results  were  noted  by 

o 

Flack  et  al.  with  outer  layer  similarity  for  Q2  and  Q4  events. 


//  8 


FIG.  8.  Contribution  to  conditionally  averaged  Reynolds  shear  stress  for  all 
surfaces  in  outer  scaling:  (a)  Q2  events,  H= 0;  (b)  Q4  events,  H- 0;  (c)  Q2 
events,  H=  2;  (d)  Q4  events,  H=  2. 


C.  Velocity  triple  products  and  higher  order  moments 

_ + 

The+ distributions  of  the  normalized  triple  products  u'3 

and  v'3  ,  as  well  as  the  normalized  axial  and  wall-normal 

_ + 

turbulent  flux  of  the  Reynolds  shear  stress,  u'2v'  and 
u'v'2  ,  are  shown  in  Fig.  9.  Andreopoulos  and  Bradshaw'7 
pointed  out  that  the  velocity  triple  products  are  very  sensitive 
to  changes  in  surface  condition.  However,  due  to  their  inher- 
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FIG.  9.  Velocity  triple  products  for  all  surfaces  in  outer  scaling:  (a)  u '3  ;  (b) 

v'3  ;  (c)  u'2v'  ;  (d)  u'v'2  .  Uncertainty  in  u'3  ±16%,  v'3  ±25%, 

_ +  _ + 

u'2v'  ±17%,  u'v'2  ±19%. 


ently  high  experimental  uncertainty,  it  is  difficult  to  use  them 
to  draw  definitive  conclusions  with  regard  to  the  extent  of 
the  roughness  influence.  That  notwithstanding,  all  of  the 
higher  order  moments  display  similarity  within  measurement 
uncertainty  for  5k  and  y>3ks.  Differences  in  rough  and 
smooth  wall  profiles  are  present  in  the  near-wall  region.  As 


indicated  in  Fig.  9(a),  u'3  for  the  rough  surfaces  is  positive 
in  the  near-wall  region,  becoming  negative  in  the  outer  layer. 
The  smooth  wall  remains  negative  across  the  entire  boundary 
layer.  A  difference  between  rough  and  smooth  walls  in  the 
distribution  u’3  is  not  surprising  considering  the  destruction 
of  the  near-wall  peak  that  is  present  in  the  profiles  of  u'2  for 
a  rough  surface.  Structural  differences,  including  the  breakup 
of  streamwise  vortices  on  rough  walls,34  change  the  mechan¬ 
ics  of  momentum  transfer  near  the  wall.  Outer  layer  agree¬ 
ment  in  u'3  between  the  profiles  for  rough  and  smooth  walls 

37 

has  also  been  observed  by  Andreopoulos  and  Bradshaw, 
Bandyopadhyay  and  Watson,  Schultz  and  Flack,  and  Flack 
et  al.s 

The  aforementioned  experimental  uncertainty  makes  it 
difficult  to  draw  a  conclusion  regarding  the  near-wall  behav¬ 
ior  of  v  .  The  smooth  wall  profiles  lie  within  the  bounds  of 
the  rough  surface  profiles  near  the  wall.  Significant  differ¬ 
ences  in  v'3  have  been  observed  for  two-dimensional 

5  38  5 

roughness.  ’  The  mesh  roughness  of  Antonia  and  Krogstad 
also  differed  from  the  smooth  wall  for  v'3  ,  however,  the 
difference  was  not  as  dramatic  as  the  transverse  rod  rough¬ 
ness,  and  the  shape  of  the  profile  was  similar  to  the  smooth 
wall  in  the  outer  layer.  Schultz  and  Flack7  noted  differences 
in  v'3  out  to  8 k  for  flow  over  a  packed  bed  of  spheres,  while 

g 

Flack  et  al.  observed  differences  in  the  triple  products  only 
out  to  5k  or  3 ks  for  milder  sandgrain  and  mesh  roughness. 

The  normalized  axial  and  wall-normal  turbulent  flux  of 
_ +  _ + 

the  Reynolds  shear  stress,  u'2v'  and  u'v'2  ,  also  demon¬ 
strate  outer  layer  similarity  as  indicated  in  Figs.  9(c)  and 
9(d).  +Antonia  and  Krogstad5  show  significant  changes  in 
u'2v'  in  the  outer  layer  for  the  flow  over  the  rods,  but  simi¬ 
lar  shaped  profiles  and  near  collapse  for  the  mesh  roughness. 
Andreopoulos  and  Bradshaw'  noted  that  the  general  shapes 
of  the  u'2v'  and  u'v'2  profiles  in  the  outer  layer  were  the 
same  for  rough  and  smooth  walls,  while  near-wall  differ¬ 
ences,  including  a  change  in  sign,  were  observed  within  10k. 
They  concluded  that  the  near-wall  changes  in  the  triple  prod¬ 
ucts  were  associated  with  strong  sweeps  near  the  surface  of  a 
rough  wall. 


IV.  CONCLUSIONS 

These  results  provide  support  for  Townsend’s  Reynolds 
number  similarity  hypothesis  for  uniform  three-dimensional 
roughness  in  flows,  where  k  and  ks  are  a  significant  fraction 
of  S.  Based  upon  these  findings,  even  comparatively  large 
roughness  elements  that  have  a  significant  effect  on  the  mean 
flow,  can  be  considered  a  small  perturbation  to  the  boundary 
layer  with  regards  to  the  outer  layer.  While  previous  results 
have  shown  outer  layer  differences  between  smooth  and 
rough  walls,  the  present  results  show  that  the  roughness  in¬ 
fluence  in  the  outer  layer  is  not  observed  for  S/k=?  19  or 
S/ks^5.5.  While  one  would  be  tempted  to  test  larger  rough¬ 
ness  elements,  obtaining  smaller  values  of  S/k  is  difficult  in 
a  boundary  layer,  which  has  been  given  sufficient  develop¬ 
ment  region  for  the  velocity-defect  profiles  to  achieve 
streamwise  self-similarity.  As  the  roughness  height  increases. 
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the  growth  rate  of  the  boundary  layer  increases  as  well,  lim¬ 
iting  the  ratio  of  8/  k,  which  can  be  generated.  In  developing 
boundary  layers  or  flows  with  a  fixed  outer  scale  (e.g.,  pipes 
and  channels),  much  smaller  values  of  S/k  may  be  obtained. 
In  these  cases,  changes  in  the  outer  layer  turbulence  may  be 
expected  to  occur  as  at  some  point  as  the  flow  must  no 
longer  retain  the  character  of  a  wall-bounded  shear  flow,  be¬ 
coming  dominated  by  bluff  body  wakes.  For  boundary  layers 
that  have  reached  self-similarity,  there  does  not  appear  to  be 
a  critical  roughness  height  that  produces  modifications  to  the 
turbulence  that  are  observed  throughout  all  or  most  of  the 
boundary  layer.  Instead,  the  outer  flow  is  only  gradually 
modified  with  increasing  roughness  height  as  the  roughness 
sublayer,  extending  5  k  or  3 ks  from  the  wall,  begins  to  occupy 
an  ever  increasing  fraction  of  the  outer  layer. 
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